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Abstract 

Our  recent  studies  of  CO  preferential  oxidation  (PrOx)  identified  systematic  differences  between  the  characteristic  curves  of  CO  conversion 
for  a  microchannel  reactor  with  thin-film  wall  catalyst  and  conventional  mini  packed-bed  lab  reactors  (m-PBR’s).  Strong  evidence  has  suggested 
that  the  reverse  water-gas- shift  (r-WGS)  side  reaction  activated  by  temperature  gradients  in  m-PBR’s  is  the  source  of  these  differences.  In  the 
present  work,  a  quasi-3D  tubular  non-isothermal  reactor  model  based  on  the  finite  difference  method  was  constructed  to  quantitatively  study 
the  effect  of  heat  transport  resistance  on  PrOx  reaction  behavior.  First,  the  kinetic  expressions  for  the  three  principal  reactions  involved  were 
formed  based  on  the  combination  of  experimental  data  and  literature  reports  and  their  parameters  were  evaluated  with  a  non-linear  regression 
method.  Based  on  the  resulting  kinetic  model  and  an  energy  balance  derived  for  PrOx,  the  finite  difference  method  was  then  adopted  for  the 
quasi-3D  model.  This  model  was  then  used  to  simulate  both  the  microreactor  and  m-PBR’s  and  to  gain  insights  into  their  different  conversion 
behavior. 

Simulation  showed  that  the  temperature  gradients  in  m-PBR’s  favor  the  reverse  water-gas- shift  (r-WGS)  reaction,  thus  causing  a  much 
narrower  range  of  permissible  operating  temperature  compared  to  the  microreactor.  Accordingly,  the  extremely  efficient  heat  removal  of  the 
microchannel/thin-film  catalyst  system  eliminates  temperature  gradients  and  efficiently  prevents  the  onset  of  the  r-WGS  reaction. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  deep  removal  of  CO  in  a  EL-rich  stream  is  a  criti¬ 
cal  step  in  PEM  fuel  cell  applications.  In  order  to  prevent 
poisoning  of  the  fuel  cell  electrodes,  the  CO  concentration 
needs  to  be  reduced  from  ~1%  to  below  10  ppm  while  con¬ 
version  of  H2  is  minimized.  This  step  is  referred  to  as  CO 
preferential  oxidation  (PrOx)  (Eqs.  (1)  and  (2)).  In  our  other 
work  [1],  we  demonstrated  that  silicon  microreactors  coated 
with  Pt/AEOs  thin-film  catalyst  can  effectively  remove  CO 
to  below  10  ppm  at  180  °C.  In  that  study,  discrepancies  of 
the  CO  conversion  temperature  dependence  were  discovered 
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between  the  reaction  results  of  the  microreactors  and  those 
reported  with  conventional  lab  reactors  in  the  literature: 


1  AH=—67  kcal/mol 

co  +  |o2 - ^  co2 

(1) 

1  A  7/ =—58  kcal/mol 

h2  +  ±o2 - h2o 

(2) 

A  H=9. 8  kcal/mol 

CO2  +  H2 - >  CO  +  H2O 

(3) 

The  studies  based  on  conventional  lab  reactors  with  particu¬ 
late  catalyst  found  that  there  existed  a  narrow  operating  tem¬ 
perature  window  for  acceptable  CO  conversion  after  light-off 
[2-4],  followed  by  declining  CO  conversion  as  temperature 
was  increased  (Fig.  1).  In  contrast,  our  studies  with  microre¬ 
actors  showed  essentially  100%  CO  conversion  between  180 
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Nomenclature 

Ai  pre-exponential  factor  in  kinetic  expressions 
Da  Damkohler’s  number  for  heat  transport 
Ei  activation  energy  for  kinetic  expressions 
(J  mol-1) 

Ftot  total  reactant  flow  (m3  s-1  at  20  °C  and 
101,330  Pa) 

A//ads  adsorption  energy  (J  mol-1) 
ki  specific  reaction  rate  for  r\ ,  r'2  and  r'3  (i=  1,2 

or  3) 

K\  K J  adsorption  equilibrium  constant  for  CO 
K3  equilibrium  constant  of  the  r-WGS  reaction 
/  reactor  axial  position,  /  =  0  for  the  entrance  of 

the  reactor  (m) 

/ 7max  the  axial  position  where  the  temperature  (at 
r-  0)  is  the  highest  (m) 

L  total  reactor  length  (m) 

Pi  partial  pressure  of  reactant  species  i  (i:  CO,  O2, 

C02,  H2,  H20)  (Pa) 

r  reactor  radial  position,  r  =  0  for  the  center  of 

the  reactor  (m) 

r’x  reaction  rate  of  CO  in  CO  oxidation,  moles  of 
reactant  per  gram  of  active  metal  of  catalyst  per 
second  (mol  kg-1  s-1) 

r2  reaction  rate  of  H2  in  H2  oxidation,  moles  of 
reactant  per  gram  of  active  metal  of  catalyst  per 
second  (mol  kg-1  s-1) 

r'3  reaction  rate  of  CO2  in  r-WGS  oxidation, 
moles  of  reactant  per  gram  of  active  metal  of 
catalyst  per  second  (mol  kg-1  s-1) 

R  reactor  radius  (m) 

Sqo2  CO2  selectivity  (defined  as  the  ratio  of  O2  re¬ 
acted  with  CO  and  total  O2  reacted) 

Cat  effective  interparticle  heat  conduction  dis¬ 
tance:  catalyst  thickness  for  the  microreactor 
or  reactor  inner  radius  for  m-PBR  (m) 

T  local  temperature  in  the  reactor  (°C) 

rw?  7b  reactor  inner- wall  temperature  (°C) 

WHSV  weight  hourly  space  velocity,  defined  as  mole 
of  total  reactant  flow  per  mole  of  precious  metal 
in  the  catalyst  per  hour  (h-1)1 
Xi  conversion  of  species  i  (/:  CO,  O2,  H2) 

Greek  letter 

Pi  reaction  order  for  rate  expression  of  H2  oxida¬ 
tion  (/:  O2,  H2,  CO) 


and  280  °C  while  there  is  only  a  slight  drop-off  (<  1%)  at 
300  °C,  suggesting  a  much  wider  operating  window  for  CO 
conversion  (Fig.  1).  An  earlier  explanation  of  this  undesired 


Non-SI  unit  used  by  convention. 


Fig.  1.  CO  conversion  vs.  temperature:  comparison  of  the  results  of  mi¬ 
croreactor  to  results  of  other  PrOx  studies  in  the  literature.  WHSV:  our 
work,  1500 h-1;  [2],  1250b-1;  [31,75b-1;  [41,200b-1. 

CO  conversion  falloff  at  high  temperatures  was  the  competi¬ 
tion  between  H2  oxidation  and  CO  oxidation  [2] .  Further  in¬ 
vestigations  found  that  the  falloff  is  caused  by  the  heat  trans¬ 
port  limitations  of  conventional  lab  reactors  [3-6]:  the  fast 
surface  chemistry  of  the  exothermic  CO  and  H2  oxidations 
leads  to  the  accumulation  of  reaction  heat  in  the  catalyst  bed; 
this  situation  then  results  in  higher  local  temperatures  (hot 
spots)  and  the  favorable  kinetics  of  the  r-WGS  reaction  (Eq. 
(3)),  which  eventually  limits  the  net  CO  conversion. 

300  °C,  suggesting  a  much  wider  operating  window  for  CO 
conversion  (Fig.  1).  An  earlier  explanation  of  this  undesired 
CO  conversion  falloff  at  high  temperatures  was  the  competi¬ 
tion  between  H2  oxidation  and  CO  oxidation  [2] .  Further  in¬ 
vestigations  found  that  the  falloff  is  caused  by  the  heat  trans¬ 
port  limitations  of  conventional  lab  reactors  [3-6]:  the  fast 
surface  chemistry  of  the  exothermic  CO  and  H2  oxidations 
leads  to  the  accumulation  of  reaction  heat  in  the  catalyst  bed; 
this  situation  then  results  in  higher  local  temperatures  (hot 
spots)  and  the  favorable  kinetics  of  the  r-WGS  reaction  (Eq. 
(3)),  which  eventually  limits  the  net  CO  conversion. 

Various  researchers  have  shown  an  awareness  of  the  im¬ 
portance  of  thermal  management  of  PrOx  reactors.  Roberts 
et  al.  [6]  studied  the  r-WGS  reaction  in  PrOx  with  an  adia¬ 
batic  monolith  reactor  and  found  that  full  O2  conversion  in 
the  adiabatic  PrOx  reactor  caused  the  downstream  tempera¬ 
ture  to  increase  to  ^300  °C  with  an  inlet  temperature  of  only 
170  °C,  thus  favoring  the  endothermic  r-WGS  reaction.  Fur¬ 
ther,  a  ID  reactor  model  was  developed  by  Choi  et  al.  [4] 
to  evaluate  overall  PrOx  performance  under  various  reactor 
heat  exchange  conditions  (adiabatic,  isothermal,  etc.).  They 
discovered  that  the  net  CO  conversion  drops  significantly  as 
the  reactor  operation  changes  from  isothermal  to  adiabatic 
condition  with  the  inlet  gas  temperature  of  200  °C.  However, 
both  articles  only  considered  heat  transfer  in  the  axial  direc¬ 
tion  ignoring  radial  heat  transfer  resistance. 

In  our  additional  PrOx  study,  we  discussed  heat  transfer 
limitations  for  both  the  microreactor  and  the  mini  packed- 
bed  lab  reactors  (m-PBR’s)  based  on  Mear’s  criterion  [1]. 
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Due  to  the  poor  thermal  conductivity  of  the  porous  catalyst 
materials,  we  concluded  that  significant  thermal  gradients 
could  build  up  in  both  axial  and  radial  directions  even  for 
the  m-PBR’s  which  are  generally  considered  to  possess  low- 
radial  transfer  resistance  due  to  their  mm-scale  diameters. 
Under  the  same  reactor  inner-wall  (boundary)  heat  exchange 
conditions,  the  microreactor  coated  with  thin-film  wall  cat¬ 
alyst  removes  reaction  heat  much  more  efficiently  than  the 
m-PBR’s  because  of  its  extremely  small  catalyst  thickness 
(~5  x  10-6  m)  and  the  resulting  short  radial  heat  conduction 
distance  compared  to  that  of  the  packed-bed  catalyst  used  in 
m-  PBR’s  (2  x  10  3  m).  The  comparatively  longer  heat  con¬ 
duction  distances  of  the  m-PBR’s  thus  cause  significant  radial 
temperature  gradients  compared  to  the  microreactor. 

In  order  to  quantitatively  examine  the  impact  of  heat  trans¬ 
port  limitations  on  the  PrOx  reactor  performance,  we  created 
a  quasi- 3D  computational  model  for  both  the  microreactor 
and  m-PBR’s.  First,  the  reaction  kinetic  expressions  based 
on  the  three  major  reactions  involved  in  PrOx  (Eqs.  (1)— (3)) 
were  formed  and  their  parameters  evaluated.  Then  a  quasi- 3D 
reactor  model  was  built  with  the  finite  difference  numerical 
method  coupling  both  the  mole  and  energy  balance  equations. 
Finally,  results  from  the  model  are  discussed  to  rationalize 
the  difference  of  PrOx  performance  between  the  microreac¬ 
tor  and  m-PBR’s. 


2.  Experimental 

2.7.  Microreactors,  catalysts  and  parallel  microreactor 
test  bed 

The  kinetic  data  for  PrOx  was  collected  with  the  sil¬ 
icon  microreactors  used  in  our  research  [7].  The  silicon 
chips  used  in  this  study  were  fabricated  with  well-known 
micromachining  processes  [8].  Photolithography  and  deep 
reactive  ion  etching  (DRIE)  by  inductively  coupled  plasma 
(ICP)  were  the  major  techniques  applied.  Anodic  bonding 
of  the  silicon  chip  to  a  piece  of  Pyrex®  glass  closed  the 
reactor,  before  or  after  the  catalyst  incorporation  in  the  mi- 
crochannel.  All  reactors  discussed  in  this  paper  have  single 
channels  with  cross-sectional  dimensions  of  5  x  10  4  m 
(width)  x  4.7  x  10-4  m  (depth)  x  4.5  x  10-2  m  (length). 

The  Pt/Al203  catalyst  with  2  wt.%  platinum  was  synthe¬ 
sized  using  a  sol-gel  technique  [9,10].  The  catalyst  precursor 
was  then  selectively  deposited  in  the  microchannel.  Multiple 
layers  were  deposited  by  repeating  the  procedures  followed 
by  calcination.  A  2  x  icr6  to  5  x  10  6  m  catalyst  thick¬ 
ness  resulted  in  the  catalyst  weight  of  roughly  0.5  x  10-3 
to  1.5  x  10-3  g. 

A  microkinetic  array  for  fast  catalyst  screening  and  pro¬ 
cess  studies  developed  in  our  lab  was  used  for  reaction  tests. 
In  this  setup,  four  reactors  are  analyzed  together  with  a  gas 
sampling  valve  multiplexing  between  reactors.  An  online 
micro-GC  (Varian  4900)  and  mass  spectrometer  (SRS  QMS- 
200)  are  shared  by  the  array  for  product  analysis.  Process 


control  and  data  acquisition  are  automated  with  a  Lab  VIEW 
program. 

Details  of  the  microreactor  fabrication,  thin-film  catalyst 
synthesis  and  the  microkinetic  array  are  found  elsewhere  [1]. 


2.2.  Kinetic  rate  expressions  for  reactions  in  PrOx 


The  majority  of  the  reported  PrOx  kinetic  studies  has  fo¬ 
cused  solely  on  the  CO  oxidation  [11-13].  However,  the  in¬ 
corporation  of  rate  expressions  of  the  coupled  H2  oxidation 
and  r-WGS  reactions  is  necessary  for  accurate  representation 
of  PrOx  reaction  behavior.  Despite  the  importance  of  evalu¬ 
ation  for  all  three  reaction  expressions,  only  few  in  the  liter¬ 
ature  have  addressed  kinetic  expressions  for  all  three  PrOx 
reactions  [4] .  Based  on  our  previous  work  [1,14]  and  the  PrOx 
kinetic  studies  in  the  literature  [4],  the  kinetic  expressions  for 
PrOx  reactions  were  formed  as  in  the  following  equations: 


r  1  = 


= 


k\PcoPo2 
(\  +  KvPco)2 


r\  =  k3 


Pcp2 

K3 


-  PcqPh2o 


K\  ~  K®  exp 


^  A/Tads 


K3  =  exp 


4577.8 


T 


(4) 

(5) 

(6) 

(7) 

(8) 
(9) 


The  majority  of  the  PrOx  studies  used  the  power  law  ex¬ 
pression  for  CO  oxidation  due  to  its  simplicity  [11-13].  This 
form  is  simplified  from  a  Langmuir-Hinshelwood  (L-H)  ex¬ 
pression  and  not  suitable  for  small  CO  concentrations  [5]. 
Therefore  we  adopted  the  full  L-H  expression  instead  of  the 
power  law  expression  for  CO  oxidation  (Eq.  (4)).  The  H2  ox¬ 
idation  was  previously  modeled  using  empirical  power  law 
rate  expressions  by  others  [4].  However,  in  the  presence  of 
CO,  the  rate-limiting  CO  desorption  strongly  inhibits  H2  and 
O2  adsorption  and  hence  H2  oxidation  in  PrOx  [1].  Thus  the 
incorporation  of  Pqo  in  the  H2  oxidation  rate  expression  is 
necessary  (Eq.  (5)).  The  kinetics  of  r-WGS  reaction  was  well 
studied  previously  [15],  in  which  an  empirical  reversible  rate 
expression  [16]  is  attractive  due  to  its  relative  simplicity  and 
its  appropriateness  in  PrOx  kinetic  studies  as  demonstrated 
previously  [4] .  With  the  rates  of  CO  oxidation,  H2  oxidation 
and  r-WGS  reaction  obtained,  the  net  rates  of  individual  re¬ 
actant  species  (i.e.,  r'CQ ,  PCq2,  /jj  ,  rf^0,  t*q  )  are  calculated 
as  in  the  following  equations: 


(10) 


~rU2  ~  rH20  —  r2  +  r3 


(ID 
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-r'0l  =  0.5(ri  +  r'2)  (12) 

A  multiple  non-linear  regression  analysis  with  the  Mar- 
quardt  method  [17]  was  selected  to  evaluate  the  independent 
variables  (A  i ,  A2 ,  A3 ,  E\ ,  E2 ,  £3 ,  P 1 ,  Pi ,  ft ,  K\ ,  A Hads)  in  the 
kinetic  expressions.  Following  the  initial  input  of  the  inde¬ 
pendent  variables,  the  dependent  variables  (Pco,  Po2,  ^co2> 
Ph2,  ^h2o)  were  calculated  with  the  4th  order  Runge-Kutta 
method  [17].  Data  fitting  of  the  calculated  dependent  vari¬ 
ables  to  experimental  data  was  then  carried  out  to  derive  cor¬ 
rected  values  of  the  independent  variables  with  the  Marquardt 
method  [17].  Afterwards,  these  corrected  values  were  fed  to 
the  kinetic  expressions  for  the  next  iteration.  The  final  values 
of  the  dependent  variables  were  reached  by  minimizing  the 
weighted  sum  of  squared  residuals  for  dependent  variables. 
The  above  algorithm  was  realized  with  MATLAB. 


/•f 


x  (q+h  1) 


R 


GJ)' 


(l.D 

k— 


\ 

s 

k 

(AO 


R  =  dr  x  q 
A  =  dl  x  p 


Ve  =  2n[(j-\)dr\  *dr*dl 


I 


(1,0 

L 


(l./H-l)  / 

— H 


Fig.  2.  The  finite  difference  grid  in  radial  (r)  and  axial  (/)  directions.  Each 
grid  point  is  denoted  as  (/',/)  with  R  as  the  reactor  radius  and  L  as  the  reactor 
length. 


2.3.  Quasi- 3D  non-isothermal  reactor  model 


The  first  part  of  the  non-isothermal  reactor  model  was  the 
construction  of  mole  balance  and  energy  balance  equations. 
The  kinetic  rate  expressions  of  different  reactant  species  in 
Eqs.  (10)— (12)  were  used  as  the  mole  balance  equations.  The 
pressure  drop  across  the  reactor  was  neglected  in  the  mod¬ 
eling  for  m-PBR’s  since  the  Ergun  equation  [18]  suggested 
acceptably  negligible  pressure  drop  for  the  flow  conditions 
and  reactor  characterizations  used  in  the  modeling.  The  first 
law  of  thermodynamics  was  used  for  the  derivation  of  the 
energy  balance  equation  (Eq.  (13))  [18].  After  dissecting  the 
heat  flow  (d Q)  into  radial  and  axial  terms,  a  partial  differ¬ 
ential  equation  (Eq.  (14))  is  obtained  as  the  energy  balance 
equation  in  the  model.  The  first  bracketed  term  is  the  heat 
removal  in  axial  and  radial  directions  and  the  second  is  the 
heat  generation  from  the  three  reactions  in  PrOx: 


d Q 

dV 


5  r 


E 


i=\  L 


H?(Tr) 


d  F 
dV 


=  0 


(13) 


l  dr  d2r 

r  dr  dr 2 


5  r- 


E 


i=l  L 


H°ATr) 


d Ff 

dV 


=  0 


(14) 


The  next  part  describes  the  quasi-3D  structural  model  of 
the  reactors.  Cylindrical  geometry  was  used  to  approximate 
both  the  microchannel  and  the  packed  bed.  Thus  the  3D  reac¬ 
tor  structure  can  be  represented  by  the  quasi- 3D  finite  differ¬ 
ence  grid  in  radial  and  axial  directions,  as  shown  in  Fig.  2.  The 
volume  of  a  differential  3D  element  is  therefore  calculated 
as  27r[(j  —  l)dr]  dr  dl  (Fig.  2). 

As  the  third  step,  derivatives  were  transformed  to  finite 
difference  forms  by  a  well-established  approach  [19,20].  The 
resulting  linear  equations  are  shown  in  the  following  equa¬ 


tions: 


cj2  1) 

0  7  07) 


1 


dr 2 

(dr)2 

cj2  n r(k-{ - 1) 

°  \fi) 

1 

_  ( 

di 2 

(di l)2 

C\rp(k+1) 

1 

= 

Ik) 


(k) 


fk) 


~  rAV\i o+i,o  Zio;i)  ^  ^j-ur 


n(k) 

07+1) 


r(k) 

0,0 


n(k) 

07-1) 


dr 


2  dr 


/ nr(k )  _  ’j-’(k)  \ 

u  0+1,0  iO-i.0; 


(15) 


(16) 


(17) 


Finally,  the  Gauss-Seidel  iteration  method  [20]  was  used 
to  solve  the  mole  balance  and  energy  balance  equations  si¬ 
multaneously.  The  sum  of  squared  residuals  was  minimized 
as  the  criterion  to  end  iteration. 


3.  Results  and  discussion 

3.1.  Qualitative  analysis  of  heat  transport  limitations  in 
reactors 

Our  additional  PrOx  studies  [1]  discussed  the  heat  trans¬ 
port  resistance  with  Mear’s  criterion,  which  suggested  that 
interparticle  heat  transport  dominates  in  both  the  microre¬ 
actor  and  m-PBR’s  compared  to  intraparticle  and  interphase 
heat  transport.  The  Damkohler  number  (Da)  was  then  used 
for  qualitative  comparison  of  heat  transport  limitations  for 
the  two  types  of  reactors,  in  which  the  heat  conduction  dis¬ 
tance  (tcat)  is  the  determining  factor.  The  /cat  of  the  thin-film 
catalyst  (5  x  10-6m)  is  orders  of  magnitude  smaller  than 
that  for  the  packed-bed  catalyst  (~2  x  10~3  m).  Due  to  the 
quadratic  dependence  of  Da  on  the  tcat.  Da  of  the  m-PBR 
with  2  x  10-3  m  radius  (2  mm  m-PBR)  is  more  than  five  or¬ 
ders  of  magnitude  of  Da  for  the  microreactor.  The  Mear’s 
criterion  then  suggested  severe  heat  transport  limitations  for 
the  m-PBR’s  for  the  highly  exothermic  PrOx  reaction. 
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3.2.  Temperature  gradients  for  the  microreactor  and  the 
2  mm  m-PBR 

In  this  section,  in  order  to  derive  the  correlations  between 
reactor  size  and  reaction  performance,  we  discuss  the  tem¬ 
perature  distribution  and  PrOx  performance  for  both  types  of 
reactor  based  on  identical  operating  conditions  (reactor  wall 
temperature  and  weight  hourly  space  velocity  (WHSV)). 

The  modeling  results  of  the  microreactor  showed  essen¬ 
tially  isothermal  temperature  distribution  in  the  thin-film  cat¬ 
alyst  even  at  the  highest  operating  temperature  (300  °C),  sup¬ 
porting  our  qualitative  analysis  with  Mear’s  criterion.  Also 
as  predicted,  the  results  for  the  2  mm  m-PBR  showed  sig¬ 
nificant  temperature  gradients  and  effect  on  PrOx  perfor¬ 
mance,  as  shown  in  Figs.  3  and  4.  Fig.  3  plots  the  conver¬ 
sions  (O2,  CO,  H2)  and  selectivity  (CO2)  along  the  reactor 
length  for  three  representative  wall  temperatures  (Tw  =  120, 
180,  220  °C).  Fig.  4  shows  the  3D  figures  of  the  temperature 
distribution  also  at  these  wall  temperatures. 

As  described  in  Section  2.3  and  Eq.  (14),  the  overall  prod¬ 
uct  of  the  reaction  rates  of  the  PrOx  reactions  (Eqs.  (l)-(3)) 
and  their  reaction  enthalpies  determine  the  local  heat  gen¬ 
eration  in  the  reactor.  Since  the  reaction  enthalpies  for  the 
exothermic  oxidations  of  CO  and  H2  are  similar  (Eqs.  (1) 
and  (2))  and  the  endothermic  r-WGS  reaction  has  reaction 
enthalpy  almost  one  order  of  magnitude  smaller  (Eq.  (3)), 
the  amount  of  local  heat  generation  is  controlled  by  the  ox¬ 
idation  reactions  and  roughly  proportional  to  the  local  O2 
reaction  activity.  This  argument  allows  the  correlation  of  the 
modeling  results  of  reaction  and  temperature  distribution  as 
follows. 

Due  to  the  low  O2  reaction  activity  at  TW=120°C 
(Fig.  3a),  the  temperature  gradient  is  negligible  in  both  axial 
and  radial  directions  (<1  °C),  as  shown  in  Fig.  4a.  At  180  °C, 
O2  reaction  activity  is  mild  up  to  /  =  2  x  10-2  m,  which  leads 
to  small  heat  accumulation  and  temperature  gradients  in 
this  region.  However,  the  sudden  increase  of  both  CO  and 
H2  oxidation  activity  between  l-  2  x  10“2  to  2.5  X  10  2  m 
causes  a  dramatic  increase  of  net  O2  reaction  activity,  lead¬ 
ing  to  temperatures  much  higher  internally  than  at  the  wall 
(Figs.  3b  and  4b).  Figs.  3c  and  4c  show  the  results  of  reac¬ 
tion  and  temperature  distribution  at  Tw  =  220  °C;  more  severe 
temperature  gradients  are  developed  close  to  the  entrance  of 
the  reactor  since  full  O2  conversion  is  reached  at  a  length  of 
only  10~3  m. 

In  order  to  gain  further  insight  into  the  effect  of  wall  tem¬ 
perature  on  hot  spot  formation,  the  axial  temperature  dis¬ 
tribution  in  the  center  (r  =  0.0m)  and  the  radial  temperature 
distribution  at  /  =  /rmax  are  plotted  in  Figs.  5  and  6,  respec¬ 
tively.  Fig.  5  shows  that  the  temperature  gradients  become 
more  dramatic  and  move  upstream  with  higher  wall  temper¬ 
atures.  With  the  wall  temperature  increases,  full  conversion 
of  O2  (with  CO  and  H2)  takes  place  in  shorter  reactor  lengths. 
Due  to  the  similar  reaction  heat  of  CO  and  H2  oxidations,  the 
total  heat  generated  by  full  O2  conversion  is  almost  constant 
even  as  selectivity  varies.  A  shorter  length  for  full  O2  con- 
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Fig.  3.  O2,  CO,  H2  conversion  and  CO2  selectivity  along  the  reactor  length 
of  2  mm  m-PBR  with  Tw  =  120,  180,  220  °C.  WHSV  =  1500  h"1 . 


version  thus  implies  a  higher  density  of  heat  accumulation 
and  consequently  larger  temperature  gradients  in  the  region 
close  to  /jmax.  Correspondingly,  the  radial  temperature  gra¬ 
dients  also  increase  with  the  wall  temperature,  as  shown  in 
Fig.  6. 
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Fig.  4.  3D  plots  of  reactor  temperature  distribution  for  the  2  mm  m-PBR 
with  rw  =  120,  180,  220  °C.  WHSV=  1500  h”1. 

3.3.  Operating  temperature  window  ofPrOx 

Fig.  7  plots  the  modeling  results  of  CO  conversion  at  dif¬ 
ferent  wall  temperatures  for  the  microreactor  and  m-PBR’s 
with  2  x  10  3  and  4  x  10  3  m  radii  (2  and  4  mm  m-PBR’s), 
all  with  the  same  WHSV  (1500  h-1)  and  isothermal  wall 
temperatures.  The  CO  conversion  curve  for  the  microreac¬ 
tor  essentially  coincides  with  the  result  for  ideal  isothermal 
operation.  On  the  contrary,  the  CO  light-off  curves  for  the 
m-PBR’s  shift  to  lower  temperatures  and  CO  conversion 
drops  significantly  at  higher  wall  temperatures.  The  light- 


Fig.  5.  Temperature  distribution  in  the  axial  direction  at  r  =  0.0m  for  the 
2  mm  m-PBR  at  different  wall  temperatures.  WHSV:  1500  h_1 . 


off  shift  to  lower  temperature  is  due  to  an  increase  in  CO 
oxidation  rate  at  local  hot  spots.  With  further  temperature 
increase,  however,  these  temperature  non-uniformities  acti¬ 
vate  the  r-WGS  reaction,  leading  to  a  drop  of  net  CO  con- 


Fig.  6.  Temperature  distribution  in  the  radial  direction  at  /  =  ljmax  for  the 
2  mm  m-PBR  at  different  wall  temperatures.  WHSV:  1500  h_1 . 
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Fig.  7.  CO  conversion  vs.  different  reactor  wall  temperatures  for  the  mi¬ 
croreactor,  2  and  4  mm  m-PBR’s.  WHSV:  1500  h_  1 . 
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Fig.  8.  The  average  reaction  rates  for  CO  oxidation,  H2  oxidation  and  the 
r-WGS  reaction  at  Tw  =  120,  180,  220  °C  for  the  microreactor,  2  and  4  mm 
m-PBR’s.  All  vertical  axes  have  the  same  scales,  with  full  scale  (1.0)  corre¬ 
sponding  to  3.9  mol/kg  s. 

version.  These  phenomena  are  seen  to  intensify  as  the  ra¬ 
dial  thermal  resistance  increases  (4  mm  versus  2  mm).  The 
inability  of  the  4  mm  m-PBR  to  reach  100%  conversion  is 
indicative  of  severe  hot  spots  even  at  relatively  low  wall  tem¬ 
peratures.  The  CO  conversion  curves  from  literature  PrOx 
studies  with  m-PBR’s  [2-4]  agree  qualitatively  with  these 
modeling  results  and  thus  can  be  understood  with  the  above 
discussion.  However,  due  to  the  lack  of  detailed  kinetic  data 
of  other  PrOx  catalysts  in  the  literature,  further  effort  to 
specifically  predict  PrOx  behavior  for  these  PrOx  reactors  is 
impractical. 

Fig.  8  then  shows  the  average  reaction  rates  of  CO  oxida¬ 
tion,  H2  oxidation  and  r-WGS  reaction  for  the  microreactor, 
the  2  and  4  mm  m-PBR’s  at  the  same  reaction  conditions 
used  in  Fig.  7,  which  were  calculated  by  integrating  and  then 
averaging  the  reaction  rates  over  the  entire  reactor  volume. 


It  clearly  indicates  that  the  m-PBR’s  has  increasing  r-WGS 
reaction  rate  with  wall  temperatures  higher  than  220  °C.  Con¬ 
sequently,  these  phenomena  are  more  significant  for  the  4  mm 
than  the  2  mm  m-PBR.  In  contrast,  the  isothermality  within 
the  microreactor  effectively  minimizes  the  extent  of  the  r- 
WGS  side  reaction  at  these  temperatures. 

Additional  results  from  our  model  also  clarified  other  fac¬ 
tors  that  influence  the  reactor  temperature  profile  and  PrOx 
reaction  performance.  For  example,  perfect  wall  insulation 
leading  to  adiabatic  operation  results  in  heat  removal  only 
through  the  axial  direction  and  as  expected,  eliminates  radial 
temperature  gradients.  Heat  accumulation  is  more  severe  than 
in  the  isothermal  wall  condition,  leading  to  an  even  higher 
local  temperature  and  adverse  PrOx  performance.  Similarly, 
a  higher  catalyst  active  metal  density  causes  large  O2  conver¬ 
sion  in  shorter  reactor  lengths  and  higher  local  temperatures. 
Thus  thermal  management  through  the  control  of  reactor  ra¬ 
dius,  heat  exchange  conditions  and  catalyst  loading  density 
is  crucial  in  the  PrOx  reactor  design. 


4.  Conclusions 

Divergence  of  the  PrOx  reaction  behavior  between  a  mi¬ 
croreactor  with  thin-film  catalyst  and  m-PBR’s  were  discov¬ 
ered  and  delineated  here.  A  quasi-3D  finite  difference  non- 
isothermal  reactor  model  was  developed  to  gain  insight  into 
the  origins  of  the  differences.  The  microreactor  has  negligible 
temperature  gradients  within  the  entire  experimental  temper¬ 
ature  range  due  to  its  small  fcat.  On  the  contrary,  due  to  its 
greater  effective  tc at,  the  m-PBR’s  develop  significant  tem¬ 
perature  gradients  in  both  radial  and  axial  directions  under  the 
conditions  considered  here  experimentally  and  through  sim¬ 
ulation.  For  wall  temperatures  before  and  during  light-off,  the 
slight  local  temperature  gradients  in  the  m-PBR’s  cause  the 
light-off  curve  to  shift  to  lower  temperatures.  At  moderately 
higher  temperatures,  severe  temperature  gradients  rapidly  de¬ 
velop  in  the  m-PBR’s  which  degrade  the  PrOx  performance 
by  activating  the  r-WGS  reaction  and  decreasing  the  net  CO 
conversion. 
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